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Background 


LGMSS model will be suspended by magnetostatic REPULSION 

- rarely used before at large gaps 

- wind tunnel systems employ attraction or combined attraction/repulsion 

Spacing between model and electromagnets is >> model scale, same order as electromagnet scale 

- electromagnet fields not affected by presence of model 

- model “sees” applied fields/gradients relatively independent of details of 

electromagnet geometry 


Problems : 

1 - llow does model behave in applied fields/gradicnts? 

2 - How to (efficiently) create desired fields/gradients? 


The Large Gap Magnetic Suspension System (LGMSS) has been described in the two previous 
presentations. The analytic approach adopted is similar to that used for many years with wind tunnel 
Magnetic Suspension and Balance Systems (MSBS), which are also large-gap systems. The motivation 
for the present study is the concern that the use of a repulsive suspension approach may present new 
problems of stability and dynamics. 



Model 


Model is cylindrical 
permanent magnet, uniformly 
magnetized in the x direction 


z 



X 


z 



I s 


Allow small translations and rotations from datum 

Use Euler angles - 0 X (roll), 0 y (pitch), 0 Z (yaw) to specifiy orientation 


The cylindrical model core (magnetic material contained within the final model envelope) is the 
originally specified configuration. Axial magnetization is the most natural choice and is the 
configuration chosen for wind tunnel MSBSs. The key point of interest is the natural behaviour of the 
model in the quasi-steady applied fields required to suspend the deadweight of the model. 
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Governing Equations 

6F = ( M . V ) B 8V ST = ( M x B ) <5V 

B is applied B from electromagnets, 

Provided M, B are relatively uniform over magnetic core then : 


F = Vol ( M . V ) B„ Illroid 
T = Vol ( M X B ccnlroid ) 


For the configuration chosen, these approximate and simplified equations are adequate and again 
correspond to traditional practice with wind tunnel MSBSs [Refs 1,2] 



Force and Torque Equations 

Expanding force and torque expressions and using small rotations : 
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With small displacements : 

B x = {n x } o 4 {b xx } o , 4 {B xy } 0 ,+ {b m } o . 4 etc 
"xx = {Bxx} 0 + {(Bxx)x} o X 4 {(»„),}„» + " 


Following some manipulation, the force and torque equations reduce to this form. Further details can 
be found in References 3,4. By way of illustration, the {Vol 1% B z } term in the T y equation is 
considered the “primary” term and arises directly from the expansion of the vector cross-product. In 
the same equation, the {Vol 0 y B x } term indicates a tendancy for the magnetization vector to align 
itself with an applied field - the “compass needle” effect. 

The effect of core translations is incorporated in the evaluation of fields at the model centroid. The 
subscript “o” implies evaluation at the datum, unt ranslated origin. 


359 




Reduction to State-Space Form 

Choose model STATE of { !iy % dy % V* V y V 5 it y * } 

- no torque about x axis - roll degree-of-freedom 

- model is initially in equilibrium (determined separately) 

Write perturbation equations in STATE-SPACE form - 

{ix} = j {ix} + j (< ij 

|specify “weighting” matrix (W), cry model mas. and inertia on leading diagonal 


W th further man, pula, ion, detailed more fully in References 3,4, the equations of motion can he 
reduced to State-Space form, where th. model state is actual,, a perturbation from equilibrium The 
equilibrium conditions, notably the electromagnet currents, me determined sepcr.tely 
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The A matrix determines the system’s dynamics. The form shown above is reasonably general. The 
field and field gradient terms are the equilibrium suspension values. Thus for any desired system 
configuration and detail design, values can be quickly evaluated and the model dynamics found from 
the eigenvalues and eigenvectors. 

The ^ matrix is filled with coefficients specifying the field and field gradient components generated by 
each electromagnet in turn. 
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Iq — ^max 


5 - Coil System 

{-0.7753, -0.2'! 17, 0.6293, 0.6293, -0.2417 } 



Coil # 

B x (Tesla) 

1 B y 

»z 

B xx (T/nn) 

B xy 


1 

0.0216 

0 

-0.0198 

0.0092 

0 


2 

0.0067 

0.0206 

-0.0198 

-0.0269 

0.0118 


3 

-0.0175 

0.0127 

-0.0198 

-0.0046 

-0.0191 


4 

-0.0175 

-0.0127 

-0.0198 

-0.0046 

0.0191 


5 

0.0067 

-0.0206 

-0.0198 

-0.0269 

-0.0118 

- calculated using OPERA (from VF/GFUN, TOSCA), using cartesian polynomial fitting 


field at grid points 








Natural Modes 



MODE 1 

14.3 rad/s 

Unstable divergence 

x, 0 y 

Axial + pitch 

MODE 2 

4.6 rad/s 

Stable oscilliatory 

x, 0 y 

Axial + pitch 

MODE 3 

12.9 rad/s 

Unstable divergence 


Yaw rotation 

MODE 4 

5.0 rad/s 

Stable oscilliatory 

z 

Vertical motion 

MODE 5 

2.5 rad/s 

Unstable divergence 

y 

Lateral translation 


The electromagnet configuration corresponds to one of the final designs emerging from the Madison 
Magnetics Incorporated design study (Reference 5). Equilibrium suspension is achieved with the current 
distribution shown, e.g. 77.53% of maximum design current in Coil number 1, the sign determined by 
the arbitrary sign convention chosen (see later Figure). The matrix can be constructed directly from 
the field and field gradient terms indicated in the Table. The coefficients of the A matrix are found by 
summation of the products of each coefficient with the relevant current fraction. Field calculations are 
carried out using an analysis and post- processing package, “OPERA”, which uses numerically 
evaluated integral expressions for field around simple conductor geometries. Once the A matrix is 
found, eigenvalues and eigenvectors are found using “PC-MATLAB”. The frequency, stability and 
shape of each mode is of interest. Mode 5 is rather benign (low frequency) and represents the model, in 
a sense, “falling off 1 the suspension electromagnets. Stability in suspension height is expected in a 
repulsive mode suspension and is shown in Mode 4. Mode 3 is the “compass needle” effect, with the 
model attempting to reverse direction so as to align the magnetization vector with the axial (B^) field. 
Mode 1 appears to represent similar behaviour in the orthogonal plane, though coupling into 
translation is exhibited. Mode 2 is an unexpected result. Unstable translation would have been expected 
by analogy with Mode 5. 




5-Coil System 


Inner radius 


0.173m 



The important dimensions of the Madison Magnetics design and the sign convention for positive 
current direction are shown here. The levitation height, measured from the model axis to the top face 
of the electromagnets, is 0.9144m (36 inches). 
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Natural Modes 



The three unstable modes are 

Mode 1 - shown top left; Mode 3 - shown center and Mode 5 - lower right 

The two stable, oscilliatory modes are 

Mode 2 - shown top right and Mode 4 - shown lower left 
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Origins of Modes ? 

Suppose only B xz applied (to support weight of model) 


MODE A 9.3 rad/s Unstable divergence x, Oy Axial + pitch 

MODE B 9.3 rad/s Stable oscilliatory x, Oy Axial + pitch 

Adding B x (typical value) 


MODE A 

14.3 rad/s 

Unstable divergence 

x > °y 

Axial 4- pitch 

MODE B 

6.9 rad/s 

Stable oscilliatory 

x, 0 y 

Axial + pitch 

MODE C 

12.9 rad/s 

Unstable divergence 

K 

Yaw rotation 


In order to better understand the origin of Modes 1 and 2, the A matrix was re-solved with 
unnecessary field and field gradient terms arbitrarily zeroed. With only the gradient required to 
generate the lifting force on the model (opposing weight), two modes are found, vaguely approximating 
Modes 1 and 2 as previously shown. If the axial field is re-applied, still holding all “second-order 
gradients” (terms of the form (B-j)^ zero, the stable oscilliatory mode is moved to a lower frequency, 
the unstable mode to a higher frequency. The “com pass- needle” mode (Mode 3) now appears. 
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Fraction of full design current 


Effect of Rotation in Azimuth (Yaw) 


Variation of coil current with yaw orientation 



Yaw orientation (degrees) 

— Coil 1 
- -Coil 2 

Coil 3 

-•-Coil 4 
Coil 5 


As the model rotates about the vertical axis, electromagnet currents are smoothly redistributed between 
electromagnets. The current variation in each electromagnet is virtually sinusoidal. It is found that 
there are no significant changes in the model’s modes of motion as the rotation proceeds. 



Six Coil System 


Try to control B x ( =>0 ). 


Find 


B 

B 


0.933 

1.0 

1.0 

0.933 

1.0 

1.0 


Currents are all large and roughly +/-/+/-/+/- for normal suspension with B x 0 


In an attempt to lower the frequency of the highest frequency unstable mode, an attempt is made to 
control the value of B x , preferably forcing it to zero. This is only feasible if an additional eletromagnet 
is added to the configuration. This is done by preserving the same individual electromagnet geometry 
and spacing between electromagnets, but locating all at a larger radius. 

It is quickly found that the resulting configuration is ineffective in controlling B x . This is due to the 
fact that the current distribution required to generate the gradient B xz is virtually identical to that 
required to generate B x , each field or field gradient component being generated independently of all 
others in both cases. 
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The cause of the problem is easily understood with a North-South pole representation. The model 
requires a N-S pole pair distributed along the x axis as shown. This inevitably generates an axial field 
in the direction opposite to the model magnetization. 

It can be noted that in wind tunnel MSBSs, this problem either docs not arise, if electromagnets are 
located symmetrically above and below the model, or is exploited, where electromagnets are only 
positioned above the model, creating a useful natural magnetizing field. 
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Assuming that the axial field cannot be sensibly controlled, attention is turned to the possibility of 
generating a lateral gradient of the “lift” field gradient, B xz ^, which has been exploited in wind tunnel 
MSBS work in the past, for generation of rolling moment. If possible, this would provide a means for 
controlling the 6th degree of freedom, presently presumed to be passively stabilized. 


It is found that the system does not work for the 6-electromagnet arrangement, due to the “roll” field 
being generated by a current distribution identical to that required for the generation of lateral force, 
fy. Additional electromagnets are therefore added, in this case requiring a reduction in size of each, 
along with an increase in the radius of their centers. The total ampere-turns in each electromagnet is 
held constant. It is found that the problem of inseparability of roll and sideforce again arises with the 
8-clectromagnet arrangement. With only 7 electromagnets, the symmetry of the arrangement is of a 
sufficiently low order to permit separation of these two fields, though not very effectively. 
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Future Work 


Incorporate modal analysis into control system simulation and design 

- may require coupled axial and pitch degrees-of-freedom 

- high frequency unstable modes place burden on power supplies and controller 

in LGMSS application 
Study second-order effects 
Study influence of eddy currents 


Simulation efforts are underway, extending work reported in References 3,4. To achieve optimum 
performance, the coupling between the axial translation and pitch degrees-of-freedom needs to be 
addressed. Existing work with wind tunnel MSBSs has dealt with similar effects by insertion of a 
“decoupling” matrix into the control loops, such that the controller can be configured as 5 (or 6) 
parallel and quasi-independent loops, each stabilizing one of the natural degrees-of-freedom. Other 
approaches are possible. 
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Conclusions - N-coil ring 


Will work for 5 degree-of-freedom control 

360 degree azimuth ( 0 Z ) range easily achieved 

Large motion capability around datum expected 

Have to control all currents in (superconducting) electromagnets 

6th degree-of-freedom can be passively stabilized if ~g present 
Alternative roll control schemes are available 

Mg present, highest frequency modes are a problem 
With {» absent, modes may not be a problem 


These conclusions are based on this study and the results of previous design studies for 5- and 6- 
elcctromagnet configurations. At present, design work for the LGMSS is focussing on an alternative 
configuration, where the axis of magnetization is vertical, parallel with the gravity vector. 




Speculation 


Change the direction of magnetization ?? 

With vertical magnetization ( parallel to ~g ), B zz 
to provide roll/pitch stability. 

Modes would be : 

Stable oscilliatory 0 X 

Stable oscilliatory Oy 

Neutral 0 Z 

Unstable divergent x 

Unstable divergent y 

??????? i 


would support weight. Can ( must ) arrange B z 

Provided B z is correct sign 
Provided B z is correct sign 
Depends on system axisymmetry 
Highly dependant 
on electromagnet 
configuration 


It appears possible to achieve a lower value of the frequency of the highest frequency unstable modes, 
by aligning the magnetization vector with the gravity vector. Further analysis is required of this and 
other possible configurations. 
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